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ABSTRACT 

This document describes the mathematical model developed to describe the 
three-dimensional motion of the Dynamic Docking Test System (DDTS) active 
table. The active table is modeled as a rigid body supported by six 
flexible hydraulic actuators which produce the commanded table motions. 
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1.0 INTRODUCTION 

The development of the three-dimensional mathematical model and computer 
program which simulates the dynamic motion of the DDTS active table in 
response to table motion commands is documented in two volumes. Volume 1 
presents the derivation of the mathematical model and Volume 2 (Reference 
1) describes the resulting computer program, "NASA Advanced Docking System 
(NADS)." 

The active table shown in Figure 1-1 is modeled as a rigid body, and each 
of the six actuators is modeled as a flexible rod with pinned ends. The 
model includes nonlinear hydraulic equations for the hydraulic actuators 
and a mathematical representation of the electronic control system for each 
actuator. Actuator position, velocity, and differential pressure across 
the hydraulic piston are used as feedback signals in the control system. 

The nomenclature used in the equations is shown in the Appendix. 



/ 



1-1 




D2-118544-1 












D2-118544-1 


2.0 COORDINATE SYSTEMS AND TRANSFOR MATIONS 

2.1 INERTIAL COORDINATES (Xj, Yj. Zj) 


The inertia! simulator coordinate system Is an orthogonal right h nd 
coordinate system whose origin is on the simulator centerlin • P la " e 

of the floor swivel joints. The y, and Zj axes form a horizontal 

plane, and the Xj axis is positive down (see Figure 2-1). 


2.2 TABLE COORDINATES (Xj, y T > z T ) 

The table coordinate system is an orthogonal, right-handed coordinate 
system whose origin is at the center of gravity of the simu ator tab 
The y T and z T axes lie in the plane of the table, and the x T 

is positive "down" (see Figure 2-1). 

2.3 ACTUATOR COORDINATES (x^ , y s .> z s .) 

Each actuator has its own coordinate system. The axis is colinear 

with the actuator centerline. The y s . axis is perpendicular to the x s . 

axis and the inertial gravity vector. The ^ axis is perpendicular to 

both x and y and is positive "up" as shown in Figure 2-1. 
s i i 


2.4 TRANSFORMATION FROM INERTIAL TO TABLE COORDINATES 

Euler angles shown in Figure 2-2 are used to transform from inertial 
coordinates to table coordinates. 

/ 
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Figure 2-1. Active Table Coordinate Systems 
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2.4 (Continued) 



Figure 2-2. Euler Angles 

The order of rotation is e, ip, <j> which corresponds to rotation about 
the y, z, x axes, respectively. The transformation from table to 
inertial coordinates is: 



Therefore: 



(2) 
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2.4 (Continued) 
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Where: 


[A] 


" Ce*Cifj 

-Cep • C0 • Sip-^S 9 • S<p 

S4> • C© • Stp+Cep *S6 “ 


Sip 

C<{> - Cip 

—Sep • C«p 

(4) 

-S0*c^ 

C(p • S0 • Sip+S<p * C6 

-Sp »S6 • Sip+Ctp *C0 _ 


C = cosine 


S = sine 




2.5 TRANSFORMATION FROM ACTUATOR TO INERTIAL COORDINATES 

The transformation from actuator coordinates to inertial coordinates uses 
the following angles: 

a, - the angle between the horizontal plane through the floor 
joint and the actuator x axis (Figure 2-3) 

P 

- the angle between the inertial Zj axis and the' projection 
of the actuator x g axis in the yj - Zj plane (Figure 2-3) 


sin e 


x i 


A i ' ' V 




cos e 


A, 


r e 2 + r 2 
s s 

y i z 


/ 


l 


Pi 


(5) 


( 6 ) 
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Figure 2-3. Actuator Transformation Angles 
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2.5 (Continued) 


sin £ 




A i 


Vr * * r * 
y i z i 


( 7 ) 


cos 6 




yjr . 2 + r 2 


(8) 


*1 


, and r are vector components of actuator length. 


where r_ , r , — - 

2 i >1 Z 1 

1 , in the inertial coordinate system, 


The transformation from actuator to inertial coordinates is then: 


*1 

*1 

z T 



(9) 


where: 


K] 


- s \ 

0 

- Ce A i 



\ 5 \ 

c % 

■ S0 A i - Se A i 


(10) 

v c \ 

A 

- C V Sfl A. 
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2.5 (Continued) 

The equivalent "Euler angle" rotations for this transformation (from 
inertial to actuator coordinates) would: 

a. Rotate -90° about yj 

b Rotate e, about the intermediate z axis 

A i 

c Rotate -e, about the intermediate y axis 
A i 

2.6 TRANSFORMATION FROM INDIVIDUAL ACTUATOR COORDINATES TO TABLE COORDINATES 

Using the previous transformations, the transformation from individual 
actuator coordinate systems to the table coordinate system becomes: 


X T 

1 

x \ 
s i 


Pr V 

y e ( 

^T 

= [ T i] 

1 



z_ \ 

Z T 


s il 


where 

W ■ [*T hi 


Li kewise: 


where : 


'■< ft! 

X T J 


y? I 

S) 1 

Z T ' 


/ 

[ T i] = [ TI i] T[A] 


(13) 


(14) 
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3.0 TARI.E MOTION COMMANDS 

Table commands are specified In the inertial coordinate system. Actuator 
commands for two types of commands will be discussed: sinusoidal position 

commands and constant velocity commands, 

3.1 SINUSOIDAL POSITION COMMANDS 

be the amplitude of commanded sinusoidal table 


Let ( 

AXj 

and | 

Ae \ 


j^I 

* 

j 


i AZ I 


A4> ) 


motion in the inertial coorainaw sy>«= - 

then obtained by adding the commanded sinusoidal motion to 

inertial position of the table. 


x T = x i + A *i sin “c 1 
c o 

y T = y, + A Yi si " “c* 

c o 

z = Z + AZ, Sin ( 0 c t 

*c o . 

e = e + A8 sin w t 

CO c 

¥ c = * 0 + ^ sin “c* (15) 

(t, = 4 + Acj) sin w t 

T c o C 

X, = AXj 0> c COS u c t 
L C 

y T = w c cos “c 1 
c 

2 = AZj « c COS u c t 

lc 

6 = A0u> COS w t 

C ^ 

= A4>u COS u> t 

w ^ 

J = A<))OJ c COS 0> c t 


t 
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3.2 CONSTANT VELOCITY COMMANDS 

For constant velocity commands (AXj, Ayj , AZj, Ad,-W, aJ ) , . the 

total inertial commands are; 



3.3 ACTUATOR COMMANDS 

The total inertial position and velocity commands are transformed to 
individual actuator commands. 

Let [aJ be the transformation from the table coordinate system to the 

inertial coordinate system (equation 4) with the Euler angles replaced by 
the commanded Euler angles. Then the commanded inertial velocities of 
the actuator/table attachment points are: 
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3.3 (Continued) 



The commanded inertial components of actuator length are 



.(19) 


Commanded actuator lengths are then 

/ 



( 20 ) 
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3.3 (Continued) 

and commanded actuator velocities are: 
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4.0 SFRVO electronics . 

The servo electronics consist of actuator position and rate command signals; 
forward loop compensation network; servo valve; and position, rate .and 
differential pressure feedback loops as shown in Figure 4-1. The forw 
loop compensation network, the valve dynamics representation, and the 
position feedback filter are designed to be optional in the computer pro- 
gram. For example, if B<1. then the forward loop compensation network 

is not included in the simulation. Similarly, the valve dynamics an 

position feedback filter are neglected if V 1 and 


u <1, respectively. 


4.1 ACTUATOR COMMANDS 
Def i ne 



«i = Kf 


'c, * V <c, 


s i 






( 22 ) 


where 1,1, and i' c are the commanded actuator length, velocity, 

c i c i i 

and acceleration, respectively. 


4.2 POSITION FEEDBACK FILTER 

The differential equation for the actuator position feedback filter is 


If w s > 1 




If 


M 



(23) 


(24) 


-4-1 



\ . 

Figure 4-1. ; Servo Electronics Block Diagram 


■v 


D2-118544- 













D2-118544-1 


4.3 DIFFERENTIAL PRESSURE FEEDBACK COMPENSATION 

Pressure ^ 

->• - f r rr 

the static differential pressure censed b, une,u.l piston areas, 
differential equations for these filters are: 


v ■ v [V (% ‘ %) ' *«i] 


a l i = V 2 


(V*h) 


(25) 


where 


is the pressure feedback gain 


4.4 FORWARD LOOP COMPENSATION NETWORK 

The forward loop compensation network consists of a lead-lag filter with 
corner frequencies a and &. 


For e>l ; 


v*' V^v 

( r L, - "f 's , 1 "l] 


(26) 


where a, > a i » x s* X s 

H i 


" , x , and x are si 


gnals from the feedback loops 


For &<!*• 


e. =0 

fe/ = 6 ^ - a Xj - K r * s - K f x s 


(27) 


where 


K, and K are the 


i M 

displacement and rate feedback gains, respectively 


4-3 



D2-118544-1 


4.5 SERVO VALVE DYNAMICS 

The dynamics of the servo valve are represented by a single-degree-of 

freedom system with a natural frequency <o v and damping rat! o c v - 


If 







(28) 


where K is the forward loop amplifier gain and Q q is the no-load 
flow through the valve. 


If 


v 1: 





(29) 
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5.0 ACTUATOR MODEL 

Each actuator is modeled as a flexible rod with pinned ends. Hydraulic 
forces are calculated using nonlinear hydraulic flow equations and unequal 
push and pull piston areas. Actuator control system electronics are 
modeled and include differential pressure, velocity, and position feedback. 
The actuator geometry and nomenclature are shown in Figure 5-1. 

5.1 ACTUATOR MASS AND INERTIA CHARACTERISTICS 
The mass moment of inertia of the piston rod is: 


where m p is the driven mass (piston and piston rod). The mass moment 

of inertia of the entire actuator assembly about the floor pivot is: 

2 


Ia = Ifl + I + m 
A. A c p p 


(V-) 


(31) 


where: 

I = Mass moment of inertia of the cylinder structure about the 
A c floor pivot (excludes driven mass, m ). 


The effective rigid lateral mass of each actuator assembly for use in the 
equations of motion is then: 


5.2 ACTUATOR FLEXIBILITY 

The dynamic bending characteristics of each actuator are calculated assuming 
that the cylinder is rigid compared to the piston rod and that the effective 
dynamic mass is lumped at the rod end seal of the cylinder. The bending 
characteristics are also assumed to be identical for each of the two bending 

planes of the actuator. 


5-1 




02-118544-1 


5.2 (Continued) 

The effective dynamic mass lumped at the cylinder end is approximated as 


% °[\ < yHv 2 ] 


(33) 


Assuming pinned joints between the cylinder and piston rod, the piston 
rod stiffness is: 



3 (EI) r l r 
1 2 ~~T 2 ”” 
rl i r2 i 


(34) 


where 1 . and l,.- are defined as follows: 
rl i r£ i 


lpl i = lp i ' lc 


1 r2 i = V ' Vl., 


The effective lateral stiffness of the actuator with a rigid cylinder is 

s H(v c -v P j)(v + ^)(vA) 


1 2 1 2 1 
‘rl. , r2 i ‘c V 


(35) 


The actuator bending frequency is then: 


= S ' \ 


(36) 


5.3 HYDRAULIC FLOW EQUATIONS 

The nonlinear hydraulic flow equations are based on the derivations presented 
In Reference 2 for a double-acting hydraulic piston. A schematic of the 
hydraulic servo valve and actuator is shown in Figure 5-2. 
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Figure 5-2. Hydraulic Servo Valve Schematic 
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5.3 (Continued) 

The flow continuity equations are: 

V 

Qi = c ( Pl - p 2 ) + h + rji 

V 

V 2 

q 2 = c p ( Pl - p 2 ) - c 6p P 2 - - i; p 2 

where: 

Ql = Q 0 - 2 K c P! 

Q 2 “ Qa + 2 K c P2 


( 37 ) 


(38) 


and 

Q = The no-load flow of the valve 

K c * Valve pressure flow coefficient 

C p = Leakage coefficient across the piston 

C = Leakage coefficient past the piston rod seal 
6 P 


The volume-stroke relationships are: 


v i = \ + A i (1 P - V 

V 2 = \ ' A 2 (l p ' 1 o ) 

= A 1 V 

^2 = " A 2 ^p 


( 39 ) 


where V and V are the hydraulic volumes at zero stroke. 
°1 °2 
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Therefore, neglecting piston rod seal leakage, the hydraulic flow equations 
for each actuator are: 

Pi-r [ Q 0 - 2K c p i ' c p (p i ' ^ ' 

1 (40) 

Pa = [' Q o ■ 2K c p 2 + C P (P 1 ■ Pz) + A 2 
5.4 ACTUATOR FORCES 

Actuator force, (f^) .re c.W« f- the dlTf«nti .1 pressure 

across the piston. In addition to the viscous d..»1»9 forces associated 
with the actuators, coulomb friction is also included. 

F p ^ A i p i - A 2 P 2 - B p 1p - C F F f (41) 


A velocity "bandwidth" for coulomb friction is used to prevent a discontinuity 
at zero velocity. 


A Force 



The coulomb friction force is: 





-C F F f 


(42) 
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5.3 (Continued) 

where C p is a coefficient which is a function of actuator velocity: 


If lipl s. v bw * then C F = 

P IV 


I f lipl < v bw ■ then C F = i®- 

bw 


(43) 
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6.0 EQUATIONS OF MOTION 

Table and actuator equations of motion are written in the body fixed 
table coordinates in the following form: 

Ixl = [M]- 1 {cl (44) 

where: |x| is a column of accelerations for each degree of freedom (six 
degrees of freedom for the table and two elastic degrees of 
freedom for each actuator) 

[M] is the 18 x 18 coupled mass matrix 

{C} is a column of generalized forces for each degree of freedom 


The mass coupling effects of the actuators due to table motions are derived 
by Lagrange's method.- The three-dimensional rigid motions of the actuators 
are completely constrained (i.e., they are dependent upon the motions of 
the table). These constraints are expressed by the velocity substitutions 
in the energy expressions. 

The final equations are much simplified when compared with the equations 
Which would result from a rigorous derivation. Due to the nonorthogonality 
between actuator and table motions, a large number of nonlinear velocity 
coupling terms results. All of these terms were assumed negligible since, 
for expected table velocities, they are quite small and their omission 
prevents the equations from becoming unwieldy. 

6.1 MASS MATRIX 

The kinetic energy of the rigid table and actuators is: 

-V 'V . 6 • 

T - *£ niy (fj • fy) + o>y • (I-y wy) + h nipl^lp^ 

6 1 (45) 

+ % E f 1 (y a 2 + * a 2 ) dm- 

i=i 0 i \ 
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6.1 (Continued) 

The first two terms in this expression are the energy of the rigid table; 
the third term is the energy of the piston rod due to stroking; and the r 

fourth term is the energy of the actuator assembly due to rotation about 
the floor pivot. 

Using the transformation from table coordinates to actuator coordinates, 
[T.] , shown in equation (14): 

V, " T <n ( ‘ " z X * "> "“i ) 

* T 1 Vi.,) 

(46) 

* T <13 ( ^ ' ‘V + 



+ T. 


12 



+ T- 


13 



a i 


■ [ T , 2 , ( V - ", V, * ", r ,., ) 


+ T, 




, 23 ( V - ", X * -x V,)] ^ 


+ Tx 


(47) 


. + V. (x s } 


= V r-^e.. *1 (x s ) 


i Pi 
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6.1 (Continued) . 

2 = T . f JL - w r + w r 1 

% L ’si v T 2 ya i y za W • 

+ T. (Vt + wT - w r i 

’32 V T 2 xa i x 2a J 

+ T* ( z T - to r + a) r — (48) 

1 33 \ T y xa i x 

+ % *i (x s } 


where: 

(j> - (x s ) is the actuator bending mode shape as a function of 

x . and y and z ' are bending velocities of the 
s e.j e.* 

actuators. 


The elastic bending modes of each actuator are assumed to be a simple 
mode shape normalized to unity at the upper end of the cylinder (x g = l c ). 

There are two identical modes for each actuator. The generalized mass 
for each mode is assumed to be lumped at the upper end of the cylinder; 
thus, the mass distribution terms can be integrated. 


e.g 


P -5 2 . 

x dm - I. = mass moment of total actuator assembly 

0 S i 


r 

about the floor pivot. 
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6.1 (Continued) 


and: 


lp 

/ 1 x s > (x s ) dm^ l c m 


V 

where m is the generalized mass of 1— actuator for each 
bending mode. 

d 3T 

Lagrange's equation requires the determination of ^ where Qj ts 

the j— generalized coordinate in the equation of motion. In this simu- 
lation: 


Ql = Xf 
^2 = y t 
^3 = Z T 

Q 4 = “ T , 

4 

Qc = “t 


Q 6 = “■ 


Q 7 = y e 3 

% = z e. 


«17 “ y e f 
^18 " z e. 


(49) 


Then: 


= * "v dr (? t • V + * t ^ ■ (I ' V] 


a Qj 


a Qj 


a Qj 


(50) 


6 3l p,- . £ pPi / 3y ai . 3Za i . \ , 

+ ■03-^1- + E J [-T-y a. + X" Z a.) dm i 

p i=19Q. P i i=l 0 ' 9 Qj 1 3 Qj 1 ' 


/ 
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6.1 (Continued) 

For j < 6: 

where C v are coefficients from equation (47). e.g. 
y ij 

Cy n = Ti 2i 

Cy i4 = rza i + Ti 23 *V a 1 
etc. 

For j = 7, 9, 11, 13, 15, 17: 



ay 


a i _ 




*. = !. 


and, for j = 8, 10, 12, 14, 16, 18 


__L = o 

3 Qj 


Likewise: 



= 0 
= 1 


for j < 6 

for j = 7, 9, 11, 13, 15, 17 
for j = 8, 10, 12, 14, 16, 18 


( 51 ) 


(52) 


(53) 


(54) 

(55) 

(56) 
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6.1 (Continued) 

Therefore, for j < 6, the last term in equation (50) becomes: 




where, from equations (47) and (48): 

V t < 2 iV t ’kVVi 

Z| *1 ° Tl 31 * a 1 + Ti 32 * a 1 + Ti 33 Z *t 


(59) 


(60) 


and, from equation (46): 


• • 


x a i = X T ' “z r ya i + “y r za i 


K. = *T + “z T xa 4 ■ “x r za 


(61) 


(62) 


z a- = Z T " “y r xa i + “x r ya i 


(63) 
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6.1 (Continued) 

Also, for j = 7, 9, 11, 13, 15, 17, the last term In equation (50) 
becomes: 


X = (y r . r- *i ^ + \ *i 2 (x s } ) dm i 

30. \ 1 P: 1 ' 



Likewise, for j = 8, 10, 12, 14, 16, 18: 



( 66 ) 


Differentiating equations (50), (58), (65), and (66) to obtain |j t ^ ^ 

J 

for j < 6 
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6.1 (Continued) 






) 


( 67 ) 


+ m 










for j = 7, 9* 11, 13, 15, 17: 



and, for j = 8, 10, 12, 14, 16, 18. 



At this point, it is convenient to redefine several quantities to simplify 
the notation. - 
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6.1 (Continued) 
also: 


Sj " T, !j 


S, ' \ 


Then equations (46), (59), and (60) become: 


\ " j^’lk ^kl 
* r i = ^ Ti 2 k ^ 


z -Y T- r 
r i 1 3k a ki 


Then, for j = 1, 2, 3: 


3Qi 


1 = y T. r + T. r 

Pi ’lk a ik 1 lk a ik 
^ r i = Sk Sk + Ti 2k fa ik 

Zr i = ^3k *" a i k + Ti 3k r “lk 


(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 
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6.1 (Continued) 

Therefore, for j = 1, 2, 3: 



The general equations are extremely complex, particularly because of all 
the centrifugal and coriolis acceleration terns. These terms can be shown 
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6.1 (Continued) 

to be small (less than 0.01 g) for the expected table velocities 
Neglecting these terms, equation (79) becomes (for j - 1. 2, 3): 



Likewise for j - 4, 5, 6: 
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6.1 (Continued) 


_ t y% - T . - r = C 
,n. ” hi ya i 1 12 xa i XZ i 


Def i ne : 


S 4 ~ ^22 rza i + Ti 23 fya i ’ CyX i 
%4 " " Ti 32 Fza i + Ti 33 ZX i 

Cy i5 ^21 *" za i ^23 xa i yy i 
Cy i6 " • Ti 21 fya i Tl 22 X3 i ^ 

Then, using equations (72), (73), and (74), equation (81) becomes: 

d_ ^ = ^row (j-3) of [ I ] t ujl + ^ X (I • S T ) j 


6 


Ik a ik 


+ ^ ( Cy ij S Ti 2k Sk + Sj S^Sk \k ) 


6 1 / .. ** \ 

£\t 


Simplifying equation (68), for j = 7, 9, 11, 13, 15, 17: 
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6.1 (Continued) 


" % S * '-1 W’a ?1 <* 


(90) 


and 


, simplifying equation (69), for j = 8, 10, 12, 14, 16, 18: 


Hf) 


1 3 

m z + m n— T. r 
*»1 e i q i 1 P 1 3k a ik 


(91) 


Expanding some 


of the summations in equations (89), (90), and (91): 


3 

V T- 

iA lk “ ik 


r. - T 


in ( Xt " fya i + Wy *** 


+ T i 12 ( *T + - r 


w z xa. 


,) 

- “x r za . ) 

+ T il3 ( Z T * “y X + “ x Vi ) 

§- T v f V v'*,) 

* t ,„ ( v* •* s - ^ s ) 


T i 23 ( Z T " “y r xa i + “x Vt) 


(92) 


(93) 


^"Sk^ik Ti 3l( 


x T - 


“z r y 3i + “y r za. ) 
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6.1 (Continued) 


* T, 3J ( y ' ' '" z r “( " “* r «i ) 

* T % ^ S * * v, ) 19,1 


Using the definitions in equations (82) through (88), these equations 
reduce to: 


k=l 


T i r a 
n lk a ik 


^i X T + + T j Z J + ‘ X X “x 

ill I i 12 I 13 1 xx i x 


+ C.... a).. + C_ u) 


(95) 


xy, “V xzi “z 


S Tl 2k ^ik 


= t t 21 x T + T i 22 y T + T i 23 Z T + c y Xi “x 


+ C oj + C ti) 

yy i y z 


(96) 


S T i r a = T i X T + T i + T i Z T + C zx. “x 
tel ’3k a ik ’31 T 1 32 1 n 33 1 zx i 


(97) 


+ C zy i “y + c z Zi “z 


The mass matrix shown in upper triangular form shown in Figure 6-1 is 
obtained by combining equations (95) through (97) and (82) through (88) 
with equations (80), (89), (90), and (91). 

6.2 GENERALIZED FORCES / 


The generalized forces are obtained by considering the work required to 
produce a unit displacement for each degree of freedom. 
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6.2 (Continued) 


fiW j = Qj {«qj) 


(98) 


where: 

0. - Generalized force 

6q. = Unit displacement 

Let F p be the net piston force along the local x $ actuator axis. The 

work done is then: 


But, since 


then: 


< w i = F Pi X s. 


(99) 


V 


■H 


( 100 ) 


6W i 


F Pi L T in Ti i 2 Ti l3^ 


T 


( 101 ) 


The generalized force for the table translational degrees of freedom are 
obtained by letting x, = 1, y T = * T = 0 ^d then y, = 1 , x, = z, = 0, 

etc. 
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6.2 (Continued) 



The displacements at the table swivel joints due to rotations of the 
are : 


AX,. 


0 

-A6 Z 

A0 y" 

I 

r 

xa. 

^1 


A9 Z 

0 

-A0 X 


| X 

AZ i 


;A6 y 

A6 X 

0 


lP »i 


For Ae x = 1, A6y = A6 Z = 0: 


Ax i 


0 

AYi 

= , 

-r 

za. 

AZ. 


r M _ 

1 

T 

ya^ 


For A 0 y = 1, " A9 X = A6 Z = 0: 


(ax. 


r za. 

ay( 

IS 

1 

0 

-r 

> 
N 
— ». 


xa i 
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6.2 (Continued) 

For &e z - 1» A 0 X = A 0 y = Os 





“1 



^i 

AZ. 

, = 

s 

1 

T 

0 


(106) 


Transforming these displacements to the servo actuator coordinate system: 



r «i 


S 

for 

Ae x = 

1 



r z a i 

-\ 3 

r -i 

for 

A3 y = 

1 

(107) 

= -T. 

r 

+ T. 

r va 

for 

A6 z = 

1 


hi 

r ya. 

h2 

xa i 






6W 

6 

1-1 

f p, 




(108) 



These generalized forces are' combined with the (S T X r T ) terms from 
equation (80), the S T X (I T • S T ) terms from equation (81), the actuator 
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6.2 (Continued) 

damping and stiffness terms and the external forces and moments to obtain 
the J C } matrix of equation (44): 
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7.0 ffll nil AT I ON OF ACTUATOR VELOCITIES ftND POSITIONS 


Actuator lengths and velocities are calculated fro, the equations of -notion 
variables for use in the servo loop feedbacks and for the determi nation of, 

actuator friction forces. 


Actuator velocities are calculated by first determining the velocities of 
the actuator/table attachment points in table coordinates: 






( 111 ) 


These velocities are then transformed to actuator coordinates to obtain 



( 112 ) 


Actuator lengths are calculated by first obtaining the components of actuator 
length in the inertial coordinate system. 


where y, and z f are the inertial coordinates of the floor swivel joints 
f i T i 




of each actuator. 
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7.0 (Continued) 


Then the actuator lengths are calculated as follows: 



( 114 ) 


These actuator lengths, and velocities are used in the feedback loops in 
the servo electronics shown in Figure 4-1. 
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Symbol 

[A] 

A2 

a. 


B P 

|C} 

C P 

e 

El.. 


EXT 


A/ 

I 


l AC 


APPENDIX 

NOMENCLATURE 

Description 

Transformation matrix from table to inertial coordinates 
."Push" and “pull" stroke working areas of actuators 

; 

Output of the pressure feedback first-order lag filter 

Output of the pressure feedback high-pass filter 

Viscous damping coefficient of actuator 

Column of generalized forces for equations of motion 
solution 

Leakage coefficient across piston seals 

Output of the forward loop compensation network 
Bending modulus of piston rod 

External force 

Coulomb friction force of actuator 

Total hydraulic and friction forces acting on pistons 

Net forces on actuator piston 

Inertia tensor of the active table 

Mass moment of inertia of entire actuator assembly 
about floor pivot 

Mass moment of inertia of cylinder (excluding the 
mass of the piston) about floor swivel joint 

Mass moment of inertia of the piston rod 
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Symbol 


! xx’ *yy ’ 


1 zz * *xy* 
X xz- ! yz 





g 

pf 

r 

rc 


1 


c 





Description 

Moment and products of inertia 


Effective lateral stiffness of the actuator 
Piston rod stiffness 
Valve pressure flow coefficient 
Displacement feedback and command gain 
Electronics and valve forward loop gain 
Pressure feedback loop gain 
Velocity feedback loop gain 
Velocity command gain 

Distance from floor swivel to center line of 
piston rod seal at end of cylinder 

Commanded actuator length 

Commanded actuator velocity 

Commanded actuator acceleration 

Retracted length (between swivel joints) of 
actuators 

Maximum stroke of actuators 
Length of piston rod 
Actuator length 
Actuator velocity 
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Symbol 


P 

1 






p r p 2 
Q j 

Q 0 

Q r Q 2 




[TI] 

t 

T 


Description 
Actuator acceleration 

Effective rigid lateral mass of actuator assembly 
Mass of piston rod and piston 

Effective bending mass lumped at rod seal of cylinder 


Table mass 

Mass matrix and mass matrix inverse 
External moments 

Moment acting about table c.g. from hydraulic and 
friction forces 

Supply pressure 

"Rush" and "pull" actuator hydraulic pressure 
jitl generalized coordinate 


No-load valve flow 


Hydraulic flow into and out of. the actuator 
Inertial vector components of actuator length 


X axis table station of actuator swivel joints with 
respect to the table c.g. 


Y, Z table coordinates of swivel joints with respect 
to the table c.g. 


Transformation matrix transforming vectors from table 
coordinates to local actuator coordinates 

Transformation from actuator to inertial coordinates 
Time ' 

Kinetic energy of the system 
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Symbol 



V z e 

y f . Zf 


a 

3 

& A 


6 

AX.j i Ayj, 
AZj, A0 » 
Aij>, A4> 

0, <t> 

e A 


Description 

Velocity bandwidth for coulomb friction 

Initial hydraulic volumes of push and pull strokes 
of fully retracted actuator 

“Push" and "pull" hydraulic volumes 
Inertial coordinates 
Actuator coordinates 
Table coordinates 

Initial inertial coordinates of table c.g. 


Bending displacements of the actuators 

Yand Z inertial coordinates of floor swivel joints 


Break frequency of first order filter 

Break frequency of first order filter 

The angle between the inertial Zj axis and the 
projection of the actuator x $ axis in the y r Zj 
plane 

Equivalent hydraulic system bulk modulus 


Total actuator command signal 

Sinusoidal amplitudes of translational commands for 
table c.q. and of table Euler angles 


Euler angles 

The angle between the yj-z.j plane and the actuator 
x g axis ' 
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Symbol 




Description 

Initial Euler angles of the table coordinate system 
with respect to the inertial system 

Actuator bending mode shape 
Damping constant for actuator bending 



Damping constant of second order filter on displacement 
feedback 


Damping constant of valve dynamics 


w 2 


Break frequencies of first order filters 


“f 

Vl’ V2 



Displacement command signal frequency 
Actuator bending frequency 

Frequency of sinusoidal .external forces and moments 

Break frequencies of pressure feedback filters 

Frequency of second order filter on displacement 
feedback 

Frequency of valve dynamics 
Table rotational rates 


/ 
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